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Q ! Abstract 

& . 

Q ' We discuss linearized gravitational perturbations of higher dimensional spacetimes. 

For algebraically special spacetimes (e.g. Myers- Perry black holes), we show that there 
exist local gauge invariant quantities linear in the metric perturbation. These are the 
higher dimensional generalizations of the 4d Newman-Penrose scalars that (in an alge- 
braically special vacuum spacetime) satisfy decoupled equations of motion. We show 
that decoupling occurs in more than four dimensions if, and only if, the spacetime admits 
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{3JT)| a null geodesic congruence with vanishing expansion, rotation and shear. Decoupling of 

electromagnetic perturbations occurs under the same conditions. Although these con- 
£NJ . ditions are not satisfied in black hole spacetimes, they are satisfied in the near-horizon 

geometry of an extreme black hole. 

o 

o : 1 Introduction 

The last decade has seen increasing interest in the study of General Relativity in d > 4 space- 
time dimensions. Black hole solutions are of particular interest [T], and the classical stability 
of such solutions is of obvious importance. It has been established that the Schwarzschild 
solution is stable against linearized gravitational perturbations for all d > 4 [2]. However, 
heuristic arguments [3] suggest that Myers-Perry black holes [4] might be unstable for suffi- 
ciently large angular momentum. 

The prediction of Ref. [3J has been confirmed by a linearized stability analysis for cer- 
tain Myers- Perry black holes with enhanced symmetry. Refs. [U E] considered the case of a 
singly-spinning MP black hole, i.e., one with a single non- vanishing angular momentum. It 
was shown that, for a critical value of the angular momentum (for fixed mass), there exists a 
stationary gravitational perturbation. It was argued that this corresponds to the threshold of 
instability, i.e., black holes with larger angular momentum are unstable. Ref. [7] considered 
the most symmetrical case, a MP black hole with equal angular momenta in an odd number 
of dimensions, and demonstrated the existence of gravitational perturbations growing expo- 
nentially with time, for sufficiently large angular momentum. Instabilities of singly-spinning 
MP black holes have also been found via nonlinear numerical evolution of a perturbed black 
hole 0. 
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In spite of this recent progress, a study of the general case seems hopelessly difficult 
owing to the complexity of the linearized equations of motion, or the difficulty of performing 
numerical evolution of solutions with a large number of free parameters. 

So far, studies of linearized perturbations of higher-dimensional black holes have exploited 
isometries of black hole spacetimes, e.g. spherical symmetry of the Schwarzschild solution, or 
enhancement of symmetry of the Myers-Perry solution when some of the angular momenta 
coincide. However, in 4 dimensions, there is an alternative approach, due to Teukolsky [5| [TP] . 
which exploits the algebraically special nature of black hole solutions. It is this approach which 
renders tractable the study of perturbations of the Kerr solution. 

The Teukolsky approach is as follows. Consider a 4d spacetime, and let (£, n, m, m) be a 
null tetrad. The Weyl tensor is encoded in the Newman- Penrose scalars \l/o, ■ ■ ■ > ^4- Now, 
consider a linearized perturbation of such a spacetime. Let denote the unperturbed value 
of a-, and let denote the perturbation. In general, there is gauge freedom correspond- 
ing to the possibility of infinitesimal coordinate transformations and infinitesimal changes 
of tetrad. However, it can be shown that is gauge invariant if (and only if) £ is a re- 
peated principal null direction of the background spacetime. Therefore, for perturbations of 
algebraically special spacetimes, there exists a local, gauge-invariant quantity, linear in the 
metric perturbation. It is natural to exploit this fact when studying perturbations of such 
spacetimes. 

In a general spacetime, the linearized equations of motion will lead to coupled equations 
for the quantities . Remarkably, in an algebraically special vacuum spacetime, Teukolsky 
showed that one can decouple these equations to obtain a single, second order, wave equation 
for . If the background is type D (for example, the Kerr metric), i.e., if £ and n both 
are repeated principal null directions, then SPq and $4 both are gauge invariant and both 
satisfy decoupled equations of motion. The solutions of these two equations are related by 
identities, and knowledge of either one is sufficient to reconstruct the metric perturbation. 

The goal of the present paper is to investigate whether, and how, these gauge invariance 
and decoupling properties extend to perturbations of higher-dimensional spacetimes. It is 
clear from the outset that, for d > 4, we should not expect to obtain a full description of 
gravitational perturbations solely in terms of a single complex scalar analagous to or 
^4 ). This is because the gravitational field in d > 4 dimensions has more physical degrees 
of freedom than a complex scalar. 

The appropriate generalization of \l/o can be identified by appealing to the classification 
of the Weyl tensor in d > 4 dimensions [TT]. This uses a basis eo = £, t\ = n, e; = m^, 
i = 2, . . . d — 1, where £ and n are null and mu\ are spacelike. The higher-dimensional analogue 
of \l/o is a (d — 2) x (d — 2) traceless symmetric matrix fl^ = C 0i0 j, where C a b c d is the Weyl 
tensor. The analogue of ^4 is another such matrix, Q'- = Cmj. These quantities transform 
as scalars under general coordinate transformations. Note that the number of independent 
components of flij (or f^ -) is the same as the number of physical degrees of freedom of the 
gravitational field. 

In d > 4 dimensions, the analogue of a principal null direction is a Weyl aligned null 
direction (WAND), and the analogue of a repeated principal null direction is a multiple WAND 
[TTj . We say that a spacetime is algebraically special if it admits a multiple WAND. It is type 
D (or O) if it admits two multiple WANDs. The MP black hole is an example of a type D 
spacetime [T2] . 

Just as for d — 4, we find that flty is invariant under infinitesimal coordinate transfer- 
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mations and infinitesimal changes of basis if (and only if) I is a multiple WAND. Q\y and 

both are gauge invariant if (and only if) I and n both are multiple WANDs, that is the 
background is type D (or 0). This gauge invariance implies that, irrespective of decoupling, 
these quantities are natural objects to consider when studying gravitational perturbations of 
higher- dimensional algebraically special solutions. 

We will study linearized gravitational perturbations of algebraically special spacetimes 
satisfying the vacuum Einstein equation (allowing for a cosmological constant). We analyze 
decoupling following the approach of Stewart and Walker |13] . This exploits the Geroch-Held- 
Penrose (GHP) formalism [T4], a very useful approach for studying spacetimes with one or 
two preferred null directions, such as algebraically special spacetimes. Ref. [15] extended this 
formalism to d > 4. We find that flj^ satisfies a decoupled equation in an algebraically special 
vacuum spacetime with d > 4 if , and only if, I is geodesic and free of expansion, rotation and 
shear. We also analyze the simpler case of a Maxwell field and find that exactly the same 
condition is required for decoupling in this case. 

A spacetime admitting a null geodesic congruence with vanishing expansion, rotation and 
shear is known as a Kundt spacetime. It was shown in Ref. [H] that any such spacetime 
is algebraically special (in vacuum). Hence our result is that decoupling occurs if, and only 
if, the spacetime is Kundt. This result can be contrasted with d — 4, for which decoupling 
requires only that I be geodesic and shearfree. These conditions are equivalent to I being a 
repeated principal null direction (by the Goldberg-Sachs theorem), i.e., they are satisfied in 
any algebraically special spacetime. 

In order for and Sl^p both to satisfy decoupled equations, £ and n both must be 
geodesic with vanishing expansion, rotation and shear. We will refer to such a spacetime as 
doubly Kundt. A doubly Kundt spacetime must be type D (or O). 

Unfortunately, black hole spacetimes are not Kundt and therefore decoupling does not 
occur in higher-dimensional black hole spacetimesrl Obviously is it disappointing that decou- 
pling does not occur for the Myers-Perry solution^ However, as we will explain at the end 
of this paper, the near-horizon geometries of known extreme vacuum black hole solutions are 
doubly Kundt solutions. Therefore our decoupled equation is ideal for studying perturbations 
of near-horizon geometries. In a companion paper [T7] we will demonstrate that, under cer- 
tain circumstances, one can predict an instability of the full extreme black hole geometry by 
using our decoupled equation to demonstrate an instability of the near-horizon geometry. 

This paper is organized as follows. In Section [2] we describe our notation, and give a 
brief summary of the necessary results regarding the higher-dimensional generalization of the 
GHP formalism from Ref. [15]. In Section [3] we investigate the existence of gauge invariant 
quantities, and show that are gauge invariant if and only if the background spacetime is 
algebraically special. We then move on to consider decoupling of perturbations. As a warm- 
up exercise, in Section H] we consider the decoupling of Maxwell perturbations, as this simpler 
example illustrates the approach that we use in the gravitational case in Section Finally, 
in Section E] we discuss the possible applications of our results. 

1 There is no contradiction with the results of Ref. [5] since that reference studies perturbations by exploiting 
the spherical symmetry of the Schwarzschild solution rather than its type D property, i.e., fiy is not used 
to describe the perturbation. The quantities that satisfy decoupled equations are non-local in the metric 
perturbation. 

2 Nevertheless, we emphasize that the quantities fiS , fiy should be useful in studies of Myers-Perry 
perturbations because of their locality and gauge invariance. 
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2 Notation 



We will make heavy use of the higher- dimensional generalization of the GHP formalism, which 
was developed in Ref. [15]. This section gives a brief summary of this formalism. 
We will use a basis eo = £, t\ = n, = mm, i = 2, . . . d — 1, where 

£ 2 = n 2 = £ ■ mm. = n ■ = 0, £ • n = 1, mm ■ = <5y. (2.1) 

Changes of basis are described by Lorentz transformations. The Lorentz group is generated 
by the following transformations: 

• Boosts (A a real function) 

£ \—> X£, n i — y X" n, mm h-> mm. (2.2) 

• Spins (Xij G SO(d — 2)) 

£ h> £, mm. i — ^ Xijim^y (2.3) 

• Null rotations about £ (zj some real functions) 

£i-> £, n i->- n + Zjm(j) — ^z 2 ^, mm !->■ mw — Zj£ (2.4) 

• Null rotations about n (z, some real functions) 

£ i — ^ £ -\- Zjm(j) — |z 2 n, n^n, mm h-> m.(j) — ^n. (2.5) 

We define 

Lab = V b la, N ab = V 6 n a , Mab = V 6 mm a , (2.6) 

and also set 

Pij — -^ij, Tj = K i — LiQ. (2-7) 

We can decompose pij into its trace p, its symmetric traceless part and its antisymmetric 
part Uij. These represent the expansion, shear and twist of the null geodesic congruence 
defined by the vector field £. /«, vanishes if, and only if, this congruence is geodesic. 

The idea of the GHP formalism is to maintain covariance with respect to boosts and spins. 
We say that an object T^...^ is a GHP scalar of spin s and boost weight b if it transforms 
as a Cartesian tensor of rank s under spins (12. 3p . and T^...^ — > A 6 !^...^ under a boost ()2.2|) . 
The quantities pij,Ti,Ki are GHP scalars with b = 1,0,2 respectively. Other quantities, e,g, 
Liq are not GHP scalars because they transform inhomogeneously under boosts or spins. 

Given any GHP scalar T, we denote by T" the object obtained by exchanging £ and n in 
its definition, so for example t[ = N i0 . This priming operation leads to a significant reduction 
in the number of equations that need to be displayed explicitly. 

The heart of the GHP formalism is a set of derivative operators that map GHP scalars to 
GHP scalars. They act on a GHP scalar T ili2 ... ia of spin s and boost weight b as: 

\>Ti 1 i 2 ...i 8 = (£-d)Ti li2 ,,, ia — bL 10 T ili2 ,, Aa + '^2Mi r oTi 1 ...i r _ 1 ki r+1 ...i a , (2.8) 

r=l 

P Ti\i2...i s = iP , '^) r Li x i 2 ...i s bLuTj^ j a -\- ^ Mi r lTii,..i r -iki T +i.~is ; (^'^) 

r=l 

^iTjih-ja = ( m (i)'d)Tjij2-ja ~bLi i T jl j 2 _j a + ^ Mj r jTj 1 ...j r _ 1 kj r+1 ...j s - (2.10) 

r=l 
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The spin and boost weights of these objects are (s, 6+1), (s, b — 1) and (s + 1, 6), respectively. 
We expand the Weyl tensor C a b c d in the frame (12. ip . and write 

Vtij = Coioj, Q,^ = Ciiij (2-11) 

^ijk = Coijk, ^ijk = Cly'fe) = CoiOfc, ^ = ClOli (2-12) 

and 

$y = Coilj, 3>yfcJ = CjjfcZ, $ = CoiOl, $^ = ^Coijj. (2.13) 

These quantities are all GHP scalars, and the notation reflects the boost weights: Q has 6 = 2, 
\l> has 6 = 1,$ has 6 = 0,^' has 6 = -1 and tt' has 6 = -2. Note the identities 

Vta = 0, Qij = ^l(ij), (2-14) 

%i = = % fe]) % i;fc] = (2.15) 

and 

— 2^ i J k 3 = ^ik = ^(ifc)> = ^[»j]> ^ = = ~~ (2.16) 

The dynamical content of GR is encapsulated in the "Newman-Penrose" equations, Bianchi 
identity and the commutators of GHP derivatives. These are written out in Appendix [A] for 
the case of an Einstein spacetime: 

R^ = Ag^. (2.17) 

We shall consider only Einstein spacetimes in this paper. 

Finally, for a Maxwell field strength F a b, we use the notation 

(pi = F 0i , F = F 01 , F ij = F ij , <p' i = F li . (2.18) 

These quantities are all GHP scalars and (pi has 6 = 1, F i and F^ have 6 = and ip[ has 
6 = -l. 

3 Gauge- invariant variables 

We are interested in linearized perturbations of spacetimes. For a quantity X, we shall 
write X = X^ + X^ where X^ ^ is the value in the background spacetime and X^ is the 
perturbation. Following Ref. [13], we look to find variables that are gauge invariant under 
both infinitesimal coordinate transformations and infinitesimal changes of basis. 

Let X be a spacetime scalar. Then, under an infinitesimal coordinate transformation with 
parameters we have X^ — > X^ + £ • dX^°\ Hence X^ is invariant under infinitesimal 
coordinate transformations if, and only if, X^ ) is constant. 

In the case of gravitational perturbations, we are interested in the scalars fly since these 
are the higher- dimensional generalization of the 4d quantity \l/o- We have 

Lemma 1 fij^ is a gauge invariant quantity if and only if £ is a multiple WAND of the 
background spacetime (or equivalently, if and only if ^fj k = = Q^- ). 

Note that this Lemma does not require any assumptions regarding the matter content of the 
spacetime. 
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Proof: First we consider infinitesimal basis transformations. Consider an infinitesimal spin 
of the form (12. 3p . If fly is non- vanishing then this will induce a change in fly^. Hence we 

must have fly = for fly to be gauge invariant. 

Next consider an infinitesimal null rotation about n. Using Ref. [T5], the change in 
is, to linear order in the infinitesimal parameters z l , 



n$»n$-2z k (*f6 j)k + *w k ) (3.i) 

+ = 0- (3-2) 

Taking the trace on j and k gives = 0. We then use = \i^(ij)k ~ ^(ik)j) to deduce 
that *g> = C 
implies that 



For invariance, we need 



that \F-°], = 0. So we conclude that invariance of fl.-J under infinitesimal basis transformations 



ng } = = o- (3.3) 

It is easy to see that these conditions are both necessary and sufficient for fly to be invariant 
under infinitesimal basis transformations. These conditions are equivalent to the statement 
that £ is a multiple WAND of the background geometry. 

Finally, since fly = 0, it follows that fly is invariant under infinitesimal coordinate 
transformations. □ 

Similarly, flA is gauge invariant if, and only if, n is a multiple WAND. Hence both 
quantities are gauge invariant if, and only if, the spacetime is type d|1 

Now consider a Maxwell field. We shall consider only a test field, i.e., we neglect gravita- 
tional backreaction and treat the Maxwell field as an infinitesimal quantity that vanishes in 
the background. It follows that all components are invariant to first order under infinitesimal 
coordinate transformations and infinitesimal basis transformations. Note that, since we are 
treating the Maxwell field as infinitesimal, and working to first order, there is no distinction 
between Maxwell theory and Maxwell theory with a Chern-Simons term. 

So far we have discussed only infinitesimal basis transformations. However, sometimes 
one might want to consider finite transformations. For example, consider a type D spacetime. 
Then i and n are fixed (up to scaling) in the background by the requirement of being multiple 
WANDs. But there is not preferred way of choosing the spatial basis vector mu\. Different 

choices are related by finite spins. flg } and Sl'V are not invariant under finite spins. Exactly 

the same issue arises in 4d, where and ^4 pick up phases under finite spins. 

Physical quantities should not care about the choice of spatial basis vectors so such quan- 
tities must be related to GHP scalars with zero spin. For example, in an asymptotically flat 
4d spacetime, the energy flux in ingoing and outgoing gravitational waves is related to the 
spin-0 GHP scalars l^o | 2 an d | 2 , respectively (for appropriate choices of £ and n, see 
[TO]). For d > 4, the analogous quantities are fly fly and fly fly . We can also define 



additional invariant quantities such as $y^fly^o 



3 Note that further gauge invariant quantities exist for higher dimensional spacetimes satisfying additional 
restrictions, see Lemma [5] later. 

4 In 4d this quantity vanishes because $^ = 5$%. 
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4 Decoupling of electromagnetic perturbations 

4.1 Main Result 

The highest boost weight components of the Maxwell field are denoted </?j. In 4d, the quan- 
tity analogous to satisfies a decoupled equation of motion in an algebraically special 
background. We shall investigate the conditions under which {pi satisfies a decoupled equa- 
tion of motion in d > 4 dimensions. The motivation for doing this is mainly that the Maxwell 
field illustrates the arguments that we shall also employ in the gravitational case, but the 
equations are considerably simpler. 

In this section, we show how, in a particular class of background Einstein spacetimes, we 
can construct decoupled 2nd order differential equations for a Maxwell test field. We show 
that this decoupling is possible if and only if the background spacetime is Kundt, that is it 
admits a geodesic null vector field that is not shearing, twisting or expanding. 

We will show that the dynamics of a Maxwell test field on the background of a Kundt 
spacetime can be described by the following equation: 

(2Kb + &A + p'b - Ar.bj + $ - fEf A ) <Pi + (-2rA' + 2r A + 2 ®l + ^tj)V 3 = 0- (4-1) 

We also show that analogous decoupled equations cannot be constructed for spacetimes that 
are not Kundt, and discuss briefly whether any alternative progress can be made. 

It is interesting to compare this to the equation of motion for a massive scalar field <fi: 

(V^V - /i 2 )0 = 0. (4.2) 

When written out in GHP form in a general background, this equation is 

(2Kb + \\ + p'b - 2tA + pb' - p 2 )4> = 0. (4.3) 

To compare this with the decoupled Maxwell equation, one must specialize to a Kundt space- 
time, for which p = 0. Note that r[ does not appear in either equation. 

4.2 Derivation of main result 

The Maxwell equations for a 2- form field strength are |15j : 



bi(pi + bF = r-ipi + pijFij - pF - Kiip't (4.4) 

2b[i(pj] - \>Fij = 2r^j] + 2Fp[ij] + 2F [i \ k p k \ j] + (4.5) 

2y Vi + b j F ji -b i F = (2p\ ij] -p , 5 ij )c Pj -2F ij r j -2Fr i + (2p m -p5 ij )<p , j (4.6) 

b {i F jk] = <p {i p' jk] + ip\ iPjk] (4.7) 



A further three equations can be obtained by priming equations fl4.4p . fl4.5l) and (14.61) . We will 
often make use of the combination &j (|4.4p — (14.51) : 

\>{Fij + SijF) = 2b[i(pj] - 5ijb k (f k - 2r'.<p i] - 2F p^ - 2F^\ k p k \^ - 

+ Sij(T' k (p k + p M F k i - pF - K k (f' k ) (4.8) 
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Now consider the combination b (l4.6p + 5, (l4.8p . This gives 

= (2b / b + 6 j 6,-)^ + 2[b,b> i -[b,6 J -](F ij + F^) + [6 i ,6,# j 

+b( - (2p' m - p'5 t] ) Vj + 2(F ij + F5ij)rj - (2 Pm - pSy)^ (4.9) 
+6i ( - PjkF jk + pF - Tjipj + Kjip'j) + bj (2T^(fj] + 2Fp [ij] + 2F\i\ k p k \j] + 2*,^) . 

This involves second derivatives of <p, as well as of the boost weight quantities Fij and 
F. However, the latter occur in the form of a commutator [b, bj}(F^ + F5ij) and can therefore 
be eliminated. Now we consider first derivatives of Maxwell components other than (p. We 
need to eliminate these from the equation if it is to decouple. 

First consider terms involving p: 

• p acts on F and F^ through the combination \>(Fij + F5ij), which we eliminate using 
equation (14. 8p . 

• Terms involving py^ are eliminated using equation f)4.6p '. 

• Terms in which p acts on pij, and p'^ are eliminated using the Newman- Penrose 
equations f lNPTj) . flNP2]) and flNP4j) ' respectively (see Appendix [£}. 

The resulting equation is very long: 

(2bb' + bjbj + p'p + pip' - 2r' j b 3 - 2r i 5 i )^ i 

+ {-p'ij\> - 2ri5j + p'jip - pijp' + [b u bj] + 2rjbi - Pj iP ')(pj 

- K^tpj - 2<pjp ki p' jk + rjPi-jp',1, - 2ip j T j T , i + 2ip i T j T , j - rjp.ip'j,, + -PjPP'ji 

+ K jV(Fij + F5ij) + PjibjF — PkibjFj k + 2p^bjF + p k jbjF ik — PjkbiFj k + PjkbjF ik 

- FbjPji - FjkbjPki + FbjPij + FijbkPjk + Fb { p - F jk bip jk + 2F\> k,i + 2F ij p'n j 
- SFpijTj - 2F^>i - AFijPjkTk - Fij^j - FjkKip'jk + F jk Kjp' ik - F jk ^ jk i - F i:j K k p' jk 

+ FijKjp' kk + FpjiTj - 3F jk pjiTk - FpTi + 2F jk pjkTi - F jk pijT k + F^prj 

+ KjbiVj + Kibj/j - Kjbjif'i + 2pkiPjk^'j + KjTjifi'i + p.kpl. jr'j 

+ f'jPikPjk - KjTif'j - KiTjip'j - Pjipif'j ~ PjkPkj^'i + 2Q ij tp' j = (4.10) 

The only terms above involving derivatives of Maxwell components other than ifi are of the 
(schematic) form htp'F, hzbip' and pbF. We need to eliminate all of these from our equations 
if we are to obtain a decoupled equation for (p^. Consider first the former two, which are 

K jP ' (F^ + F5ij) + 2^-6^;.] + Kib^ = 2 Kj p' (F i:j + F<y + . . . (4.11) 

where we have used (14. 8p ' to eliminate the nbtp' terms in favour of np'F and some other terms 
not involving derivatives. 

Now, the Maxwell equations cannot be used to eliminate the terms of the form up' F 
without re-introducing 1-derivative terms of the form ubip'. Hence the only way in which 
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the k\>'F terms can be eliminated is if = 0, hence the vector field I must be geodesic for 
decoupling to be possible. We assume henceforth that this is the case. 
Now examine the pbF terms above. These are: 



Pji&jF — pki&jFjk + 2pij5jF + pkj&jFik — PjAFjk + PjkbjFi 



ik 



(4.12) 



To achieve decoupling, we need to eliminate these terms from the equation without introducing 
any 1-derivative terms (unless the derivative acts on <p). It is convenient to decompose biFj^ 
into parts that transform irreducibly under SO(d — 2): 

2 



F ijk 



d- 3 



r8i\j&\lFi\k], 



(4.13) 



where Fijk is traceless and can be decomposed further into objects transforming irreducibly 



according to the Young tableaux 



and P. The quantity b{F, transforms in the same way 
as 6jFji, i.e. as a vector (□) under SO(d — 2). The latter can be eliminated in favour of the 
former using equation (I4.6p . which gives bjFji = biF + . . ., where the ellipsis denotes terms 
in which derivatives act only on if. The contribution of the "vector" terms to (14. 12p is then 

(Pji + {d- 3)py - pSij) bjF (4.14) 



d-3 



We can substitute our decomposition of bF into the Maxwell equations. There are no Maxwell 
equations that can be used to eliminate biF without reintroducing new derivative terms of 
the form p\xp'. Hence the only way in which the Maxwell equation will decouple is if the 
expression in brackets in (14. 14[) vanishes. The symmetric and antisymmetric parts of the 
resulting equation give 

a ij = = {d-A)u} ij> (4.15) 

where a and u are the shear and rotation of i respectively (i.e. they are the symmetric 
tracefree and antisymmetric parts of pij). Hence a necessary condition for decoupling is that 
I be shearfree and, for d > 4, rotation free (and hence hypersurface orthogonal since £ is 



geodesic). We now assume d > 4, so we set <jj 



Pij 



P 



d 



Ui 



henceforth, and therefore have 



(4.16) 



A spacetime admitting a null geodesic congruence with vanishing rotation and shear is called 
a Robinson- Trautman spacetime if p ^ and a Kundt spacetime if p = 0. It was shown in 
Refs. [lEl [TS] that an Einstein spacetime of either of these types is algebraically special, with 
the vector field i aligned with the congruence being a multiple WAND. Therefore we can take 
Qij = ^ijk = 0. Note that (INP3I) now implies bip = 0. 

It is now guaranteed that we can use equation (14.61) to eliminate "vector terms" of the 
form bF from (14.101) . Upon doing so, we find that the terms involving Tijk all drop out. The 
commutators [b, b'] and [5j,5j] can be used to tidy up the equation, giving 



= [2Kb + + P'\> + ^IpK - Arjbj] ^ + 2(r,A - rfij)^ 



+ 



2p 



3^- - *a - fz^p'm 



d-2 
d-4 



2d 3 . i . 
~~d~^T ) S - 



d-2 



P 



Tj ( Fij - F5ij) + j^ti 



(4.17) 



9 



The only term involving ip^ is the final one, so for (p[ to decouple we need (d — 4)p = 0. This 
also ensures that the terms involving Fy and F drop out of the equation. Hence decoupling 
requires p = (since d > 4), which implies p^ = 0, so £ must be free of expansion as well as 
shear and rotation. That is, the spacetime must be Kundt. The equation reduces to 

(Pi + (2Tj6i - 2Tibj + 3$ 4j - = 0. (4.18) 

which is equivalent to (14. ip . 

To summarize, for d > 4, ipi satisfies a second-order decoupled equation if, and only if, £ 
is geodesic with vanishing expansion, rotation and shear, i.e., if, and only if, the spacetime is 
Kundt. 

Note the presence of factors of (d — 4) in several of our equations above. When d = 4, 
it is not necessary for the rotation of £ to vanish in equation ( I4.15p . or for the expansion 
p to vanish in equation (14.17p . Indeed, in 4d, all that is required is that £ be geodesic 
and shearfree, which is equivalent (by the Goldberg-Sachs theorem) to the spacetime being 
algebraically special. 

It is clear that cp^ will satisfy a second-order decoupled equation (the prime of the above 
equation) if, and only if, n is geodesic with vanishing expansion, rotation and shear. Hence 
<Pi and p\ both will satisfy second order decoupled equations if, and only if, = k' { = pij = 
p[j = 0. We call such a spacetime "doubly Kundt" . 

We have considered the conditions for decoupling of a Maxwell 2-form field. It would be 
interesting to consider the more general case of a p-form field with p > 2. The equations of 
motion for such a field are given in GHP form in Ref. [To] . 



2b'b + bjbj + p'p - Arjbj + $ 



2d -3 
d-1 



A 



4.3 The Schwarzschild solution 



Consider the special case of the higher-dimensional Schwarzschild solution, which is not 
Kundt. This solution has pjj = j^^ij anc ^ r i = (a consequence of spherical symmetry). 
The latter implies that the terms in Fy and F drop out of equation (I4.17P leaving us with an 
equation of the form 

id - 4) . 



(C<p)i + 



rP>; = 0, 



(4.19) 



(rf-2) 2 ' 

where £ is a second order differential operator. The second term remains an obstruction to 
decoupling. For the Schwarzschild solution, the two multiple WANDs have identical properties 
so we can take the prime of the equation to obtain 



Jd—W p ^ 



o. 



and hence 



a ( j2 c v 



(<*-4) a 
(d-2 ' 



2 

ZP'~<Pi 



0. 



(4.20) 



(4.21) 



So in fact <pi does satisfy a decoupled equation but it is fourth order in derivatives. Note that 
we had to make use of several special properties of the Schwarzschild solution to obtain this 
result. It would be interesting to investigate more generally the circumstances under which 
one can obtain a decoupled equation of higher order for <pj. 
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5 Decoupling of gravitational perturbations 



5.1 Introduction and main result 

We now move on to gravitational perturbations. In Lemma [T] we found a set of gauge in- 
variant quantities fly under the assumption that t was a multiple WAND of the background 
spacetime, so we shall consider gravitational perturbations of an algebraically special Einstein 
spacetime, for which we can take I to be a multiple WAND. Hence fly and ^y^ vanish in the 
background, so we can treat them as first order quantities: fly = fly , \I/y-fc = ^yjL Therefore 
we shall not bother including a superscript W on fl or ^ below. 

The final result will be similar to that of the electromagnetic perturbations; we will find 
that we can only achieve decoupling when the spacetime is Kundt, that is when it admits 
a non-expanding, non-shearing, non-twisting null geodesic congruence. We will show that 
gravitational perturbations of such a Kundt spacetime are described by 

(2Kb + hh + p'p - Qr k b k + 4$ - gA) fly 

+ 4 (t*8«| - T (M b k + + 4$£, k ) Q k{j) + 2<$> ikjl Q kl = 0, (5.1) 

where all quantities except A are evaluated in the background geometry (e.g. $ denotes $^°- ) 
etc.) 

In a doubly Kundt spacetime, fly also will satisfy a decoupled equation, which is given 
by priming the above equation. 

5.2 Derivation of main result 

We follow as closely as possible the approach of Ref. [TJ]. We start by obtaining an equation 
in which second derivatives act only on fly. Consider the equations 

= -6 fe fly - 5 jk biQu + bjVL lk - b(%5 jk + V ijk ) + 5 jk ($u - 2$ a - 

+ 2(*[jfc|5ij + ^i6[ k \i + ^i[ k \i + V[k\ii)pi\j] + {^ur'iSjk - QikTj + fly-r^) (5.2) 

and 

= -2b'fly + b k (%6 jk + V ijk ) + (-flyp' + 2fl ife pj jfc] ) - 4(* (i 5 i)fc + * ( y )fc )r fe 

+ ®jkpik ~ ®kjpik + ®ikpjk ~ ®kipkj + 2®ikpkj ~ ®ijp + ®ikjipki + ®pij- (5.3) 

Equation (I5.2p is obtained by taking various linear combinations and contractions of the 
Bianchi equations (1B1I) (in Appendix |Aj, while equation ( 15. 3 j) is constructed from the sym- 
metric part of (1B2p and a contraction of (1B3|) . These equations are exact no decomposition 
into background and perturbation has been performed at this stage. 

Now we consider the linear combination 5t.f l5.2p + b( )5.3p . This contains second derivatives 
acting on fly and on \I/yfc. However, the point of taking this particular combination is that 
the second derivatives of ^y^ occur in the combination — [b, b k ](^ij k + ^iSj k ) and therefore 
can be eliminated in favour of terms involving one or zero derivatives of \&yfc using the formula 
flC2p for the commutator [b, b k ] ■ 
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We can also symmetrize the entire equation on ij without losing any useful information, 
as the antisymmetric terms do not contain any second derivatives of Q. This reduces the 
equation to 

= -(2b'b + B fc 6 fc )ny - 2[b, b']Oy - [B (ih fi fc ]n %) + WwPh) - HT^'k) 
+ [b, b k ]I* k - MT^n) - 2b {il (T* lkPkl ) + 2b l (T^ lk p m ) - 26,(T* fc p«) 

+KO«) - fiiC^li**) (5-4) 



where 



7j fc = + (5.6) 

Note that these quantities satisfy the following relations: 

Tijki — lfi\3\k)i = ^i*'|fc|0' = ^ an< ^ ^yfcj = ~~ ~~ 2)$ 4 s fc + (5.8) 

7^ = T* )fc , 7^ = and T* = \<tt> 3 (5.9) 

Tijki — T(ij) k i, = —ttijShi Tft kl = and T^ fcfe = — (cf — 2)0^. (5.10) 

In this notation, the parts of ( 15. 2 p and (15. 3p symmetric on ij become 

\> T ?jk - biTijik — -Tijik T ! + ^T{i\kiPi\3) ~ ^ T ?jiPik - ^T^ m p m i5 k \j) - T ikjl Ki (5.11) 

and 

- b k T? jk + 2\>'Q ik = T? kjlPjl - 4T* k r k + T? jkl p' kl . (5.12) 
Next we perform the following steps: 

1. Use the commutator ( 1C2[) of Appendix lAl to eliminate the terms [b, b k ]T^ k from (15. 4p 
(note that this introduces a new kind of term, of the schematic form 

2. Expand out the brackets using the Leibniz rule for GHP derivatives. 

3. Eliminate the term \>T^ k using equation (15. lip . 

4. Use the NP equations (INPip . (INP2p and (INP4I) ' from Appendix [A] to eliminate terms 
in which \> acts on pij, n and p\ - respectively. 

5. Take a linear combination of the Bianchi equations (IB2IIB3|IB4[) to get an equation 

\> T ikjl = ^' T ijkl + §{i\^l\j)k - §l^{ij)k - £(i[*tf|i)J&TO*m + SkfiyVj) 

+(-26 { * (< | + 6(i|*i)«J*|i) + (-25 fe %| + 6(i|* fc )^|i) + • • • , (5.13) 

where the ellipsis indicates terms that involve no derivatives. Use this to eliminate bT^ 
from (l5\4l) . 
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6. Use a combination of ( 1B5|) and (1B7I) to show that 




3\>'T* k + 37^,7} + . . . 



(5.14) 



where the ellipsis denotes first order terms not involving any derivatives. Use this to 
eliminate b{T® kl from ( 15. 4p . 

The resulting equation is very long so we shall not write it out in full. It has the schematic 



(b'b + 6 ■ 6 + [b, b'] + [6, 6] + p'b + p\>' + rb + r'b + rr' + pp' + 

+ k\>'V + pd^ + (tk + t'k + p 2 )$ + (rep)*' + (rp + r'p + up + b 'k + 6p)* = (5.15) 



Here, we neglect terms that are of quadratic order or higher when we decompose quantities 
into a "background" piece and a perturbation. Recall that Q and \l/ are first order quantities. 
Note that the only terms containing derivatives of Weyl components other than f2y are of the 
schematic form k,\>'^/ and p5\l/. For decoupling to occur, these must vanish for any possible 
perturbation. We shall now examine the circumstances under which we can eliminate these 



If k\ ^ then there is nothing we can do to eliminate these terms. The only Bianchi equation 
containing b'\l/ is (1B5|) . and using this again would reintroduce the 1-derivative terms that 
we have eliminated above. Hence the only way for these terms to drop out is for Ki to vanish 
in the background. Hence = is a necessary condition for decoupling. Henceforth we 
assume Ki is a first-order quantity, in which case the above terms become second order terms 
and can be neglected. 

Recall that = is equivalent to the statement that £ is geodesic in the background. 
This places no further restrictions on the spacetime, as it was shown in Ref. [19] that any 
Einstein spacetime admitting a multiple WAND also admits a geodesic multiple WAND. 

Having set k\ = 0, the only remaining terms involving derivatives of Weyl components 
other than Q are of the form pS^. The detailed form of these terms is: 



Ap{i\ h b k ^\j) + p k i [26i^f(ij) k + &(i\^i\j)k - &(i^j)ki + 2b k ^(ij)i - b(i\^ k \j)i] 

+ p ki i\ [-djtfyjifc - 8jtf/0) fc - b^k + 26 fc *y)] (5.17) 



For decoupling we need to eliminate these terms in favour of terms in which derivatives act 
only on Q. 

Certain combinations of terms of the form 6^/ can be eliminated using Bianchi equations. 
In order to understand precisely what kinds of terms can be so eliminated, we can decompose 
bi^/jki into parts that transform irreducibly under SO(d — 2). If we do the same for the 
Bianchi equations at our disposal (or combinations of them such as (15. 3p ) then we will see 
which irreducible parts of 8\1/ can be eliminated from the above equation. 

Decomposing into tracefree and trace parts gives, for d > 4: 



form 



terms. 



The detailed form of the k\>'^/ terms is 



(5.16) 



bi^jki = Vijki + 2Si[ k \Wj\i] + 5ijX k i + 25j[ k \Yi\^ + 25i[k\6j\i]Z, 



(5.18) 
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where Viju is traceless and satisfies Vi\j k \\ — V%i{M) — 0- The other terms are given by 

W M = \x i3 = 1 (-{d - 3)6 fe % ]fc + b^ j} ) , (5.19) 

% = (36**[«]* - (d - Wi]) , (5-20) 

= (rf-2)(d-4) " 3)6fe *fe> + ft (*^) " 6 ***^) ' ( 5 - 21 ) 

Y{ij) = (d-2)(d-A) ( ft **&> " (d " 3)6 ^ + 8 ***^) ' ( 5 ' 22 ) 

-5 fc ^.. (5.23) 



(d-2)(d-3) 

Note that Wry) and Y^-) are traceless and Xu^ = 0. 

The traceless part Vij k i can be decomposed further into parts that transform irreducibly 
under SO(d — 2). The relevant irreducible representations correspond to Young tableaux 
with 4 boxes. However, it turns out that we will not need to discuss these. As well as these 

quantities, we have two independent quantities transforming as B, namely W[ij] and , two 
quantities transforming as □□, namely Wuj) and Yujy and a singlet Z. 

Consider first the singlet Z, i.e., the trace 6 k ^ k . The contribution of this to equation 

dsniD is 

4{d - 3)a tJ Z = -^a^k, (5.24) 

where the shear cr^- is the traceless symmetric part of pij. In order to achieve decoupling, we 
would need to add to (15.17P a combination of Bianchi components containing a singlet term 
that cancelled this, and did not introduce any 1-derivative terms (e.g. p$ terms) that we 
have already eliminated. However, there is no such combination. For example, the singlet 
drops out of equation (15. 3p . Therefore, the only way to eliminate the singlet term from our 
equation, as required for decoupling, is to set = 0, i.e., in the background geometry, the 
shear of the multiple WAND i must vanish. Henceforth we assume that this is the case. 

Next consider the traceless symmetric tensors that arise in the above decomposition of 
bi^/jki, i.e., W(ij) and Y^y The contribution of these to ( I5.17P is: 

- 5p(W (ij) + %)) + \{d - 10) (W (ik) uj jk + W {jk) u ik ) - \{d - 2) (Y (ik) uj jk + Y {jk) u ik ) , (5.25) 

where u {j = p^y 

Now consider the Bianchi equations. The only combination of equations involving Wuj\ 
and Y(jj) that does not introduce any 1-derivative terms that we have already eliminated is 
f)5.12p . which gives an expression for 

6 k T* k = -(d-2)(W (ij) + Y (ij) ). (5.26) 

We can use this to eliminate, say, Y^ from (15.25p . i.e., we have F(y) = — W(jj) + . . ., where 
the ellipsis denotes terms in which derivatives act only on Q. ( I5.25P then reduces to 

2(d - 4) (W (ik) u jk + W w u ik ) + .... (5.27) 
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Since we have no independent equation that will allow us to eliminate W^j), we conclude 
that in order for the 5\1/ terms to decouple we must have Uij = in the background, i.e., the 
multiple WAND i must be free of both shear and rotation. Note the factor of d — 4: for d = 4, 
vanishing rotation is not necessary for decoupling^ 

Having set af^ 1 = ojf) = 0, we find that that the 1-derivative terms (I5.17P reduce to 



( 6 ^> fc + b ^ = d^2 bkT ^ k (5 - 28) 

These terms can be eliminated from ( 15.15)) with equation ( I5.12|l . 

In the resulting equation, we now use (lNP3j) to argue that bip is a first order quantity. 
It appears only when multiplied by so such terms are second order and can be dropped. 
The only Weyl components that are now acted on by derivatives are fly, and the equation 
has been reduced to the schematic form 

(Kb + 6 • 6 + [b, b'] + [6, 6] + p'b + p\>' + rS + r'b + pp + tt' + $)fl + p 2 $ + rp^ = (5.29) 

At this point, we can also simplify the form of the terms involving Q, by using the commutators 
flCl|IC3j) from Appendix |A] to eliminate the terms of the form [b, b']fl and [5, b]Q respectively, 
in favour of terms that involve at most first derivatives of Q. 

The terms of the form <3>p are simplified by noting that equation evaluated (15 .3p in the 
background geometry implies that 



p (0 H 0) = ^p (0)$( %- ( 5 - 3 °) 

Equation (I5.29P now reduces to something sufficiently simple to write out explicitly: 



(2Kb + \h + p'b + ^fpK + -hpp' - Qr k b k + 4$ - gA) a, 

+ (4r fc 6 (4 | - 4r (4 |6 fc + ^p{p' m - p' (i{k ) + 4<% + 16$^ fc ) Q m + 2<S> ikjl tt kl 

+ J^2 ~ + 2pTk " *^>* + A^**) = °- ( 5 - 31 ) 

This equation is the analogue of equation (I4.17P for the Maxwell field. Note that (I5.30p implies 
that p($fj — jZ2§5ij) is a first order quantity. To achieve decoupling we have to eliminate 
the terms not involving fly, i.e., those on the final line of this equation. For d = 4, this is 
automatic since the particular combination of $ terms appearing in this equation vanishes 
identically (i.e. = |<Mij if d — 4), as does the particular combination of ^ terms. For 
d > 4, the only way of eliminating the <3> terms above is to set p^ = 0, i.e., take p to be first 
order. All terms on the final line above are then higher order and can be neglected. 
Hence we see that, for d > 4, decoupling requires that 

«f = pf = 0, (5.32) 

i.e., the multiple WAND must be geodesic and free of expansion, rotation and shear. In other 
words, the spacetime must be Kundt. This is a necessary condition for decoupling; it is also 



5 For d= 4, &ij k = — 2<5ju*&fci , so the irreducible parts of bi^jki are just the trace, tracefree symmetric and 
antisymmetric parts of b^j. Considering the trace gives = as for d > 4. The tracefree symmetric part 
can be eliminated with (|5.3p . The antisymmetric part simply drops out of (|5.17p , using the fact that all 2 x 2 
antisymmetric matrices commute. 
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sufficient since we now have an equation in which the only perturbed Weyl components that 
appear are fiy. 

The resulting decoupled equation is: 

(2p'p + b k b k + p'b - Qr k b k + 4$ - ^A) % 

+ 4 (r fc 6 (i | - r (l ,6 fc + $| ]fe + 4$^) Q m + 2<S> ikjl n kl = 0. (5.33) 

We remind the reader that Qij is a first order quantity, so quantities multiplying Q (e.g. $, 
r) must be evaluated in the background geometry. 

5.3 Comment on the expanding case 

Just as we did for Maxwell perturbations, it is interesting to consider what happens if i is 
geodesic with vanishing rotation and shear, but non-vanishing expansion (i.e. the spacetime 
is Robinson- Trautman). Under these circumstances, we have equation (I5.3ip . a perturbation 
equation for a gauge invariant quantity f^-. However, it contains two terms that obstruct the 
decoupling of the equation. It is interesting to ask how these terms are consistent with gauge 
invariance. The answer is supplied by: 

Lemma 2 Let i be an expanding, non-twisting, non-shearing geodesic multiple WAND for 
an Einstein spacetime of dimension d > 4. Then 

(5-34) 

is a gauge invariant quantity. If 'r/ 7^ ; then 

and r< 0) *gj[ (5.35) 

also are gauge invariant quantities. 

The Schwarzschild black hole in arbitrary dimension is an example of a spacetime admitting 
such a multiple WAND (although in this case, rf^ = 0). In four dimensions, (I5.34p vanishes 
identically in all spacetimes, while the quantities (I5.35j) are not gauge invariant. 

Proof: From equation (I5.30p we have 

^ = (5-36) 

in any such spacetime. Hence we see immediately that (I5.34p is invariant under infinitesimal 
coordinate transformations, and also under infinitesimal spins. Furthermore, Ref. [TS] showed 
that all such spacetimes are of algebraic Type D so we can choose our basis so that all Weyl 
tensor components with non-zero boost weight vanish. Under an infinitesimal null rotation 
about £, equation (2.32, [15]) implies that, to first order in z iy 

$| ^ $| + z^j) - z k * (ij)k , (5.37) 

but \1/ is a first order quantity and hence 3>f- and $ = both are invariant in a Type D 
background. An identical argument applies to null rotations about n, and hence (I5.34p is a 
gauge invariant quantity. 
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For an algebraically special spacetime, $?ijk and both vanish in the background, and 
so, to first order, they are invariant under infinitesimal spins and infinitesimal coordinate 
transformations. They are also invariant under infinitesimal null rotations about £, as these 
can only introduce terms involving Qij which also vanishes in the background. We now 
consider the effect of an infinitesimal null rotation about n. Taking the prime of (2.34. [15]) 
implies that, to linear order, 

*W H. % jk + 3^%**] + z&$ k (5.38) 

and 

*i ^ *i - fci*^. ( 5 - 39 ) 

where we have used (I5.36p . We will show that the quantities ( I5.35P are invariant under this 
transformation if ^ 0. 

Take a double trace of the Bianchi equation (IB7P for the background spacetimes. This 
implies that (d - 4)5 fc $ (0) = 0, and hence, for d > 4, 5 fc $ (0) = 0. The trace of (lB5l) gives 

5 j$ (o) = grf (5.40) 

and hence = if 7^ 0. From (15.361) we then have = 0. Putting these results 

back into (B5) implies that = 0. Inserting these results into (I5.38|l5.39p implies that, 

although ^ijk is not invariant under infinitesimal null rotations about n, both r^ijk and 
are invariant, and hence both of these are new gauge invariant quantities, provided that d > 4 
and t; (0) ^ 0. □ 



6 Discussion 

We have shown that, for linearized gravitational perturbations of an algebraically special 
spacetime, there exist local quantities f^j , linear in the perturbation, that are invariant under 
infinitesimal coordinate transformations and infinitesimal changes of basis. For perturbations 
of a type D background, e.g. a Myers- Perry black hole, both and f^- are gauge invariant. 
We emphasize that, irrespective of decoupling, the locality and gauge invariance of these 
quantities should make them useful in studies of gravitational perturbations. 

We have shown that satisfies a decoupled equation of motion in d > 4 dimensions if, 
and only if, the background is Kundt, i.e., admits a null geodesic congruence with vanishing 
expansion, rotation and shear. Therefore the decoupling property which is satisfied in the 
Kerr spacetime does not extend to the Myers-Perry spacetimes dimensions in an obvious way. 

When decoupling does occur, an important question is whether a solution of the decoupled 
equation uniquely characterizes the gravitational perturbation. If one has two solutions with 
the same then do they describe the same metric perturbation? This is equivalent to 
the question of whether there exist non-trivial linearized gravitational perturbations with 

ng» = o. 

For perturbations of a Kerr black hole, this problem was addressed in Ref. [20J, where 
it was shown that a well-behaved solution with ^E'q 1 '' = must also have = 0. A null 
rotation about n can then be used to set = and a null rotation about I can be used 
to set ^3 = 0. It follows that the perturbation must preserve the type D condition to 
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linear order. Since type D solutions are specified by just a few constants [2TJ, it is natural to 
expect that there will be only a finite number of solutions satisfying these conditions. Ref. 
[2"U] showed that the only well-behaved solutions correspond simply to perturbations in the 
mass or angular momentum of the Kerr solution. 

For d > 4, even if one can show that f^j = implies that f^- = then it is no longer true 

that one can use null rotations to set \E^], = ^ ^ = 0- This is because a null rotation about 

n contains fewer parameters than the number of independent components of (whereas 
for d = 4 both have 2 degrees of freedom). So for d > 4 it seems likely that the perturbations 
overlooked by our decoupled equation are more general than perturbations preserving the 
type D condition. Nevertheless, since has the same number of degrees of freedom as 
the gravitational field, it seems reasonable to expect that our decoupled equation of motion 
captures "nearly all" of the information about linearized metric perturbations, i.e., that only 
a few special solutions are missed. 

Sometimes it might not be enough to know the solution for \ one might need to know 
the metric perturbation explicitly. Ref. [22] gives a systematic procedure for constructing 
solutions of the linearized Einstein equation (in a certain gauge), given the existence of a 
decoupled equation of motion for a quantity linear in the metric perturbation. It seems likely 
that this procedure can be applied in the present case to generate solutions of the higher- 
dimensional linearized Einstein equation whenever fljj' satisfies a decoupled equation, i.e., in 
a Kundt background. 

Finally, we discuss an application of our decoupled equation. Any extreme black hole 
admits a near- horizon geometry [23]. It turns out that all known extreme vacuum black hole 
solutions have near-horizon geometries which take the form of a fibration over AdS2 [24] : 

ds 2 = L(y) 2 (-R 2 dT 2 + ^pj + g u {y) (d^ 1 - k 1 RdT) {d<p J - k J RdT) + g MN (y)dy M dy N , 

(6.1) 

where d/dcj) 1 , I = 1, . . . ,n are rotational Killing vector fields of the black hole and k 1 are 
constants. The metric in the first set of round brackets is the metric of AdS2- The metric 
depends non-trivially only on the d — n — 2 coordinates y M . A calculation reveals that the 
vector fields dual to — dT ± dR/R 2 are tangent to affinely parameterized null geodesies with 
vanishing expansion, rotation and shear. Hence the above spacetime is doubly Kundt so our 
decoupled equations can be used to study gravitational perturbations of it. This is under 
investigation [T7] . 
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A Useful GHP equations 

Many of the results in Sections H] and [5] were based on the Newman- Penrose equations, Bianchi 
equations and derivative commutators for the higher- dimensional Geroch-Held-Penrose for- 



18 



malism defined in Ref. [15]. Here we reproduce these equations for convenience, in the case 
of an Einstein spacetime satisfying R^ u = Ag^ u . 

A.l Newman-Penrose equations 

Ippij—bjKi = -pikPkj - KiTj - TiKj - Qij, (NP1) 

Ki-P'^ = fHj(-Tj+Ti)-*i, (NP2) 
2&\j\Pi\k] = 2 T iP{jk] + lKiP'\jk] ~ ^ijk, (NP3) 
p'pij - Sj-Ti = -TiTj - KiKj - p ik p' kj - $y - jZT[8ij. (NP4) 

Another four equations can be obtained by taking the prime ' of these four (i.e. by exchanging 
the vectors I and n). 

A. 2 Bianchi equations 

Boost weight +2: 



-2(^ {jl 5u + $>i8 m + V m + Vy\u)pi\k] + 2%^ ]; (Bl) 

Boost weight +1: 

- b% - bj% + p'% = -(*^<fc-*; i fc)K Jfe + ($ ifc + 2^ + $<J i fc)p J fe i 

+(^ ljfc - *i<J ifc )^ - 2(* (< 5j- )fc + V<i j)k )T k - fi ife p' fcj , (B2) 

-\>® ijkl + 28^^ = -2*j i | fcl K| i] - 2^ mK]l] 

+4$J.p [fc q - 2$[ fc | i p i | i] + 2$^^^] + 2$ ii[fc | m p m |i] 
-2*[i|fezr' j] - 2%^] - 2fi i[fc |p^] + 2Q j[k p' m , (B3) 

— Syi^ilW] = 2®[jk\Pi\l] ~ 2®i[jPkl] + ®irn{jk\Pm\l] ~ 2&i\jP k l], (B4) 

Boost weight 0: 

- = 2(*' b .|<Ji, - *{,-|«)/?i| fc] + (2$i b -(J fc]1 - 2^$^ - $ %& )r ; 

+2(tt f %, - % W )p' m + 2%4], (B5) 

-26[<^ fc] = 2*{ < p ifc] + %,p l | fc] -2* [i ^ fc] -*, [ij -| / o(| fc] , (B6) 

-8[fc|^y|Jm] = -*i[fcJ|Pj|m] + *j'[W|Pi|m] ~ ^[k\ijP\lm] 

-^i[kl\Pj\ m ] + ^j[kl\P'i\m] - ^[k\ijP\lm]- (B7) 



Another five equations are obtained by applying the prime operator to equations (IBip -f lB5|) 
above. 
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A. 3 Commutators of derivatives 

Acting on a GHP scalar of boost weight b and spin s, commutators of GHP derivatives can 
be simplified by: 

[b, V '!):..,„ = {-Tj + 7^)8jT il „. < , + b ( ,y,j + K j ^ j + <I>) 7 

s 

+ E ('••< '•" - < K i + < T l- T ^'i + ^t) T n..,..,, (CI) 
r=l 

[b, 6<]T fcl „. fc , = + r/b + Pji&j)T kl ... ks + b (-r-pji + Kjfo + %) T kl ___ ks 

s 

+ E ( K *rPli ~ Pkri T ! + T krPli ~ P'kri K l ~ *«*r) T kl...l...k 3 , (C2) 



r=l 



[6,, » i ]T fcl ... fc . = (2p fo] b' + 24-jb) T fcl ... fcs + 6 (2^,4, + 2$£) T kl ... ks 

s 

+ (2pfc,.[<[^|j] + 2 Pk r [i\Pi\j} + $ *iM + jzih\kAj]i) T kl ...i...k s - (C3) 

r=l 

The result for [b',5j] can be obtained from ( 1C2j) / . 
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